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Abstract 

The discrete model of the real symmetric one-matrix ensemble is 
analyzed with a cubic interaction. The partition function is found to 
satisfy a recursion relation that solves the model. The double scaling- 
limit of the recursion relation leads to a Miura transformation relating 
the contributions to the free energy coming from oriented and unori- 
ented random surfaces. This transformation is the same kind as found 
with a cuartic interaction. 
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1. Introduction 

Matrix models have been used to study the formulation of non-perturba- 
tive two-dimensional Euclidean quantum gravity (pure or coupled with con- 
formal matter with c < 1) following the seminal work of references Jl|, J|, 
[0]. These models may also be viewed as models of random surfaces or zero- 
dimensional string theory. The partition function of the Hermitian matrix 
ensemble is identified as the theory of oriented surfaces. This connection is 
established by identifying the dual to the Feynman diagrams of the matrix 
model with random discretizations of two-dimensional surfaces of arbitrary 
genus. The technique of the double-scaling limit allows us to recover the 
continuum limit of the theory [|IJ , , [|J . 

Many other one-matrix models have been solved following the tools devel- 
oped for the Hermitian ensemble fl, ||, M. In this paper we shall consider 



the real symmetric one- matrix model [0], 0, [[|, ||T0[ introduced to describe 
theories with unoriented random surfaces. In these treatments, the potential 
in the matrix model was for simplicity taken to be even. When this is the 
case, the orthogonal polynomials suitable to solve the model happen to have 
a definite parity, thereby providing a computational simplification similar 
to the one appearing in the standard Hermitian model |TT]. In this letter 
we shall address the problem of a non-even potential which we solve for the 
discrete model with a cubic interaction in sections 2,3,4. 

This study is motivated by the equivalent problem taking place in the 
Hermitian ensemble where it is found fl2fl , fl3|| that the continuum limit of 
the model with a cubic potential leads to a string susceptibility satisfying 
the Painleve I equation, as it happens with a cuartic potential. In section 
5 the double-scaling limit of a real symmetric matrix ensemble with a cubic 
term is computed and compared with the cuartic potential case of 0, ||. 

Our starting point is the partition function for the real symmetric 2Nx2N 
matrices M T = M, 

Z 2N = = [ J] [ dM l\ e-^ triym) (1) 

0<i<?'<2iV-l 
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with the matrix potential given by 



V(M) = — M 2 + — M 3 (2) 
2 3 

The identification of ([!]) as a theory of unoriented random surfaces comes 
from the perturbation expansion of the interaction term in ([]]), (||) around 
the gaussian ensemble = 0. The propagator is 

< M{M n m > = g^(S im 5^ + S^J (3) 

The first term in (^) represents the "usual" propagator in the double 
line notation of the Feynman diagrams, while the second one represents a 
"twisted" propagator. The latter contribution amounts to the loss of ori- 
entability of the triangulated surface associated to the dual of a Feynman 
graphQ. Each connected vacuum diagram is weighed by a factor of N x , 
where x is the Euler characteristic of the surface, which now can take both 
even and odd values. Thus, all kinds of non-orient able surfaces contribute to 
the sum of the random surfaces described by (|3|) j/J, |J, 



Z 2 n — N 2 Z sphere + N X Z RP 2 + N°(Z torus + Z K i ein bottle) + • • • (4) 

2. The Method of Skew Orthogonal Polynomials 

The parttion funtion of non-orientable pure gravity is computed with the 
help of an appropiate set of polynomials in a similar way to what happens 
in the case of the Hermitian matrix model. Thus, when we integrate over 
"angles", i.e., over the real orthogonal group with the action M = OXO^ 1 = 
OXO T in (Tip, it is possible to express the partition function in terms of the 



eigenvalues X = diag{x\, . . . , x 2n ) of the matrix M. The result is [14| 



2N 

Z 2N = Kx J n<NA(x)|e-^-i^ (5) 

i=l 



lr To be more precise, in the continuum limit it is also needed that both contributions 
in (||) have the same positive weight 0. 



2 



where V(x) is given by (0). The main new feature is the appearence of the 
Vandermonde determinant A(x) = Yii<j{xi — Xj) in absolute value and raised 
to the power of one, rather than the usual factor A 2 (x) showing up in the 
Hermitian matrix model. This is a major difference between both models as 
far as the computation of (|5]) is concerned. From the symmetry properties 
of (|J) we may write \TA 



2N 

Z 2N = K(2N)\ { ... / IJ e -§E^(*<) x [det(x*r% )j=1) ..., 2N 

J R(-oo,x 1 ,...,x 2 n,oo) J i = i 

(6) 

here the region of integration R is — oo < X\ < . . . < x 2 n < oo. Note that 
the absolute value in (^|) is no longer required. As usual, it is convenient 
to introduce monic polynomials Ri(x) in order to rewrite the Vandermonde 
determinant as 



A(x) = \det{x) 1 )] i j=i,... 1 2N = [deti2i_i(a;j-)]ij=i,...,2jv (7) 

We want to solve the partition function Z 2 n in terms of quantities related 
to these monic polynomials. To do that, we have to chose the Ri(x) as skew 
orthogonal polynomials with respect to the following anti-symmetric scalar 
product fllCf 



/OO /"OO a 
\ dxdye{x - y)e2l v ^)+v(y)] 
-oo J —oo 

= r m Zij (8) 
where the only non-zero values of j are 

Z2i,2i+1 — ~ z 2i+l,2i = 1 (9) 

and 
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Using Mehta's method [13] of integrating over alternate variables, we 
carry out the integration over the odd variables Xi,x%, . . . ,x 2 n-i first and 
then over the remaining even variables. The final result is 



N-l 

Z 2N = K{2N)\ J] r k (11) 

k=0 

So far the problem is formally solved. In order to make the solution ex- 
plicit, the strategy consists of relating the set of skew orthogonal polynomials 
Ri with the well-known set of orthogonal polynomials solving the Hermitian 
one-matrix model. Thus, let us introduce orthogonal polynomials C, with 
respect to the follwing symmetric scalar product 

/oo 
dxe-WMCiWCjix) 
-oo 

= hAj (12) 
These polynomials are also monic and satisfy a two-step recursion relation 

xCt = C i+ i + OiCi + piCi-x (13) 

The coefficient cij takes into acount that we are dealing with a non-even 
potential (0). In addition, it is possible to derive equations for p i} Oi from 
the reltions < Cj, C i >c= and < Cj_i, C i >c= ^i-i- This yields to 

<C 1 ,PV'C 1>C =0 (14) 

<C i ^(5V'C i > c =ih^ l (15) 

Now let us express the set of polynomials Ri as a linear combination of 
the set Ci. 

Ri(x) = Ci(x) + a M _iC;_i(x) + ai A _ 2 Ci t i- 2 (x) + a iii _ 3 Ci ii _ 3 (x) (16) 

In the next section we shall see that the combination (0) ends at Cj_3 
because the order of the potential V(x) is 3. 
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Neither Ri nor C, have a definite parity anymore due to the cubic term 
in V(x). Then, it is no longer possible to split the relation flT6|) into two, one 
for the set Rn and another for i?2i-i, as occurred when dealing with even 
potentials fllP ], 0. This amounts to an extra complication in solving fll6|) . 



At the end of the computation of the coeffients cijj, we shall make such a 
splitting to notice the appearence of mixed terms between both sets. 

To proceed further, we have to relate the two scalar products introduced 
in (||) and fll^D - Given two functions /, g it is easy to prove, upon integration 
by parts, the following essential relation 

< f,g >c= -\ < pVf,g > R + < f',g> R (17) 

This equation will give us enough relations to find out a^j, provided we know 
the action of V' (x) and ^ on the set of orthogonal polynomials Cj. In fact, 
computing < Ci_j,Ri >c, j = 0,1,2,3 in fll6|) with the aid of flTTD, it is 
possible to set up the following relations for the unknown variables a^j 

= < Ci—j,R{ >c 

1 

For the action of V (x) on Cj, using (Q) and flT3|), we have 



- < PV'C^, Ri > R + < Cij, Ri > R j = 0, 1, 2, 3 Vi. (18) 



V\x)Ci(x) = g 3 C i+2 + [g 3 (<Ti+i + <r%) + 92]C i+1 

+ [9-i{pi+l +&i + Pi)+ 92CTi]Ci 

+ [93Pi{°i + Vi-x) + 92Pi]C i - 1 + g-zPiPi-iCi-2 

l=+2 

= g v i>i+l C i+l {x) (19) 

l=-2 

It is possible to simplify this expression using the discrete string equations 
(|14D and (^). In fact, the following expressions for the string equations are 
readily obtained 

= i = PlgsPifa + ffi-i) + 92pi] (20) 
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Vi,i = = g 3 (pi+i + o] + + g 2 o { (21) 

From ( |19D and fl2~0|) we can also express the coefficient f^i+i in a more 
tractable way 

= = o ( 22 ) 

Pi+l PPi+1 

On the other hand, the action of d/dx on Cj is given by 

dCj _ ^4 < Cj, Cj >c ^ 
dx j =Q < Cj, Cj >c 

which is obtained by partial integration on < Cj, >c with j < i. As V(x) 
has degree 3, this expression turns out simply to be 

dC 

= Pv ht _ 2 C t _ 2 + f3v i>i ^C^ 1 (24) 

dx 



3. Recursion Relation for the Real Symmetric One-Matrix Model 

Due to the lengthy calculations involved for establishing and solving the 
equations (|ll), it is convenient to outline the main steps needed to achieve 
the final result as follows, 

i) Choose one of the four equations in (fl"8| ) (j = 0, 1, 2, 3). 

ii) Express V' (x)Cj_j and C i _ J - in terms of the orthogonal polynomials 
d. To do this, use ©-(H) and (H- 

iii) Express the Cj polynomials of step ii) in terms of skew orthogonal 
polynomials. To do this, it is required to invert the relations fll6l). 

iv) Use the skew-orthogonality relations (H) and (^j) to solve the right 
hand side of (|T|). 

Moreover, to carry out the above program we shall make repeatedly use 
of the following essential observations, 

a) As Ri = Ci + order (Cj_i), then Cj = Ri + order (Ri_i). 

b) From the skew-orthogonality relations (§),@ we have 
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< Rlii R21+I >R — r i — — < Rfii+l, Rli > R (25) 

and the rest of scalar products are zero. 

c) The polynomial R 2 i is skew orthogonal to any polynomial of degree 
less than or equal to 2i. 

The polynomial R 2 i+i is skew orthogonal to any polynomial of degree less 
than 2i. Simbolically, 

R 2 i skew_L Ri,R 2i -i,... 

R 2i+ i skew_L Rm-u ■ ■ ■ (26) 

To be consistent, we need to show that the relations fllfiD actually stop 
at Ci_ 3 . Namely, let us assume that a term, say, Cj_ 4 is present. Then, 
the term in the r.h.s. of (|l8|) with the highest Cj is < Cj_4 + 2, Ri >r- This 
follows from ( 19 ) and (p3|). Now using observations a) and c) we easily get 



< Ci-2, Ri >r= 0. 

We shall first work out the equations fll8|) in the order j = 3,2, 1. The 
expressions for a iti _ 3 , a iti _ 2 , a iti _i thus obtained will be replaced in equation 
( |IED for j = 0. This will turn out to be the desired recursion relation that 
solves the discrete model. 

Equation ( fl^j j = 3. 

O i)i _ 3 /i i _3 = — — 1>j-3,t-l < Ri-l,Ri >R (27) 

It is convenient to split the indices into even and odd. If % = 2k, then as 

< R 2 k-i,R 2 k >r= 0, one deduces 

a 2k ,2k- 3 = (28) 
If i — 2k + 1, as < i?2fc, -R2fe+i >r= Tu and ^—3^—1 = #3, one deduces 

a2fe+l,2fe-2^2fe-2 = —^93 r k (29) 

Equation ( |7^ ) j = 2. 
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It is convenient to start with 



0-1,1-2^-2 — —^i-2,i < C i} Ri >R — < Ci-ljRi >r (30) 

After some algebra we get 

di,i-2hi-2 = —^[Vi-2,iOi,i-l + Ui-2,i-l] < Ri-liRi >R (31) 



Ui = 2k then 



d2k,2k-2 = (32) 

If % = 2k + 1 then, 

a 2fc+l,2fc-1^2fc-l = — ^[g^k+l^k + "5 ] Tk (33) 

^ PP2k 

and a2fc+i,2fc is to be determined. 
Equation j — 1. 

Applying again the procedure previously described, we obtain 
2 

— a^j-x/ii-i = — fi-i^+i < Ri + i, Ri >r 

— [^i-X,i+l( a i+l,i a i,i-l ~ a i+l,i-l) + ^i-l^i-i] < -Ri >R (34) 

Where we have kept in mind that = due to the string equation (|2l|). 

Notice that in (B3) appears aj+i^-i, which is of the same type described by 
the previous case fl31~|) . Replacing the value of a i+ i^i in fl34|) we shall obtain 
an equation for a^-i solely. However, the analysis is simplifyed if we first 
make the splitting of indices into even and odd. 
If i = 2k, then < Ri >r = and 

2 

-Tja 2 k,2k-ih2k-i = g^k (35) 
If % — 2k + 1, then < -Rj+i, -Rj >/?= and 

2 2A; + 1 

-^0'2k+l,2kh2k — —93 a 2k+l,2k[02k+2,2k+l + ^ J r £i (36) 

P PP2k+l 
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This equation admits two solutions. We shall take the simplest one, 



02fc+l,2fc — (37) 

The case a 2 k+i,2k ^ is considered in Appendix I. Now we can replace (|37j) 
into (^) to obtain the value of a 2 k+i,2k-i, 

02fc+l,2fe-l = 7 (38) 

P2k ft 2fc-l 

Now we have the solution of all the coefficients relating both sets of poly- 
nomials in (|i~6l). With this solution, it is convenient to split the relation (|16|) 
into polynomials with even and odd indices to see the picture that emerges 
from the previous calculations. Therefore, 

R2k = C 2 k + 0,2k,2k-lC 2 k-l (39) 
R-2k+l = C 2 k+1 + Cl2k+l,2k-lC 2 k-l + 0-2k+l,2k~2C 2 k-2 (40) 

where the coefficients are given by (^5|), (|37|) and fl40|) . In this way, we observe 
that there are odd contributions to R 2 k as well as even ones to R 2 k+i ■ This is 
the new effect of dealing with a non-even potential, causing the polynomials 
R2k) -R2/C+1 not to have a definite parity anymore. 

Now we are in position to set up the recursion relation between the rel- 
evant quantities r^, hk that ultimately solves our model. This is achieved 
with equation (|T8| ) for j = 0. The main steps of this lengthy calculation is 
given in Appendix II with the following result, 

h 2 k = [ — j- £ 3 T + ] r k (41) 

4 ll 2 k+l <zp2k+l 

The equation ([18]) for j = are actually two relations, depending on 
whether the index i is even or odd, as was done in the previous cases. This 
means that the system of equations is constrained, for there is one more 
equation that unknown variables. It is worth noticing that either equation 
for j = leads to the same solution (see Appendix II), thereby providing 
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a proof of the consistency of our solution. In addition, if we set = we 
should end up with the solution of the gaussian model |10| , 0, as it happens 



h 2k = ^r k (42) 

for pi = i//3 in the gaussian case. 

4. The Partition Function of the Model 

So far, we have related the calculation of the partition function for the 
real symmetric one-matrix model to the solution of the Hermitian ensemble. 
Once the norms of the orthogonal polynomials Cj are determined from the 
discrete string equations (^(]), (0), then the quantities r k are in principle 



computed from (|4l|) yielding the solution of the partition function (|TT|) . In 
fact, the recursion relation (|4l"|) obtained for r k amounts to a recursion relation 
between partition functions at different N. From ( |TTD it is possible to rewrite 
r k in the form 

Tk = (2k + 2)(2k + l) X ~2k~ ( } 

Upon substitution of ( f4"5D in pT| ) we readily obtain the following recursion 
relation for the partition function, 

^giZ 2k +4 - \(2k + A)(2k + 3)(2k + l)h 2k Z 2k+2 
(2k + 4V 

+ (2kjT h2k+lh2kZ2k = (44) 

When the continuum limit of the model is to be taken, we shall see that 
it is convenient to introduce the following quantities relted to the "norms" 
r k by 

W k = (45) 

"2fc 

and the partition function (|TlD is determined by the quantities W k and h 2k 

as 

N-l 

Z 2N = Kx (2N)\f3-( N -V ft W k h 2k (46) 

fc=0 
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Hence, the free energy takes the form 



JV-l N-l 

T 2N = In K N + J2 In h 2k + £ W 2k (47) 

k=0 k=0 

According to 0,0, the second term in T 2 n represents the contribution of 
the orientable surfaces to the free energy, for it enterely involves the constants 
related to the Hermitian matrix model which describes only oriented surfaces. 
On the other hand, the constants W k amounts to the contribution of non- 
orientable surfaces to the free energy. This is made more plausible by noticing 
that the odd powers of (1/iV) in T 2 n come precisely from W k 0. 

5. The Continuum Limit 



We now want to take the continuum limit of our solution (41 ) to describe 
the critical point of non-orientable random surfaces so, we need to borrow 
the results concerning the double-scaling limit for the Hermitian model with 
a cubic interaction. 

The scaling limit of the string equations fl20"|),([2TD goes as follows. First, 



we take the planar limit (N -> oo, with f -> 1) of ©,(121]) at i = 2N. 



Therefore, it is convenient to introduce the functions 



V x {p,&) = 2p + g 3 pa = l (48) 



V 2 (p,a) = 2a + g 3 (2p + a 2 ) = (49) 

where we have normalized g 2 = 4 without loss of generality. Then, the critical 
values of g 3c and p c are defined by 



dVi(p, a) 



(50) 



dp 

subject to 

d ^f^ =0 (51) 
dp c 

The analysis of these conditions leads to the following critical values that, 
for future convenience, we recast as 
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p\ = \ and dlcPc = ^ (52) 

Second, we insert the values (|52|) in the discrete string equations (|3)>(i3) 
at i — 2N and consider the following scaling ansatze || |7| [|T2]jf| 

P2Af+fe = Pc[l - ^ 2 /(^)] (53) 

o"2Af+fc = -5 2 s(t fc ) (54) 

where 

8 = /T 1/5 -> (55) 

t fc = (2/3 - 2JV - fc)<5 (56) 

The doble-scaling limit is achieved when 5^0 and -j| — > 1 while t = to 
is held fixed. Following reference jn§Q, it is appropiate to introduce two 



auxiliary scaling functions f±(tk) defined in terms of /?(£&) and cr(tk) by 

U(t k ) = f(t k )±s(t k ) (57) 

When dealing with non-even potentials, the key point to obtain the dif- 
ferential equations for /+ and /_ in the limit 5 — * is to express /(£&) and 
s(tk) again as a perturbative series in the S parameter^ 

f(t) = f (t)+fi(t)8 + f 2 (t)5 2 + --- (58) 

with 

/o(*) = \p-{t) (59) 

2 Let us notice that this is the same ansatze as is |l2| but in a different notation. 
3 I am very grateful to S.Dalley for explaining this point to me and drawing my attention 
to ref.Jl^] during the course of this work. 
4 We refer the reader to ref. 11131 for details. 
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and 



s(t) = s (t) + si(t)<5 + s 2 (t)5 2 + ■■■ (60) 

As far as the leading behaviour of the partition function is concerned, we 
only need to know about the differential equations satisfied the first terms in 



the expressions (pq) , (|60|) . These are found to satisfy [jT^] |0| (up to trivial 



rescaling) the Painlee I equation 



fo ~ \fo = t (61) 



suplemented by the equation 



P+\ 5=0 = fo + s = (62) 

This completes the analysis of the Hermitian model with a cubic potential. 
To begin with the continuum limit of the real symmetric model, we rewrite 
the recursion relation at i = N in the following form 

^ W N+1 W N p2N+2p2N+l - 2N 2 p 1 W N + p 2N +l = (63) 

where Wjy is given by (||5|). This is a quadratic equation in the unknown 
variable W. The spherical limit of fl5"BD is achieved by taking N, (3 — ► oo with 



N/f3 finite and assuming Wn — > W c , pn — > p c - Tuning g 3 and p to their 
critical values (p2]) previously found, we observe that ( p5| ) has a unique root 
VFc given by 

= 2p c (64) 
Let us notice though that p c can take two possible values according to 



(52|). So, the critical behaviour of the Hermitian case leads also to critical 
behaviour in the symmetric model. 

The equation is similar to the ones appearing in the analysis of the 



real symmetric ensemble with a cuartic potential , [[J and in the simplectic 
ensemble || . The technical difference is that our equation is quadratic in the 
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unknown variable W while those in the forementioned references are cubic. 
Therefore, we shall use the same scaling ansatze for W as in J5|, namely 

W N+k = W c e 6 ^ (65) 

where 5 = (3~ 1 ^, and W c , are given by (|64]), ( |56"1) respectively. The novelty 
of working with a non-even potential is that u)(tk), with t = tcb has to be 
expanded in a series of 5 

uj{t) = u (t) + + w 2 (t)5 2 + ■ • ■ (66) 

in agreement with the procedure developed in fl58|) for the functions p^. As 
far as the leading behaviour of the partition function is concerned, we are 
only interested in the differential equation obeyed by the term uj q in (|66|) . The 
way to proceed is to insert the ansatze ([53]) and fl6"5|) in the recursion relation 
flB3"D and to expand in 5 with the aid of Then, we solve for the u's in 
|66| order by order until we obtain a differential equation that determines ujq. 
This turns out to be, 

3/ = uol - 2uj' (67) 

We can thought of fl6~7|) as an inhomogeneous ordinary differential equa- 
tion yielding the solution for the unorientable contribution u to the free en- 
ergy once the orientable contribution fo is known after solving the Painleve 
I equation flBTD - This is the same kind of differential equationF] that appears 
in the analysis of the real symmetric ensemble carried out with a cuartic 
potential §, @. To be precise, the general solution found out with a cuartic 
interaction also includes in ( |6?D a term of the type ce-^i a; °( s ) ds ; where c is a 
constant of integration. Our solution corresponds to the choice c = 0. This 
c-term corresponds to different non-perturbative contributions to u> but, as 
far as the finite genus contributions is concerned, there is no difference of 
information among the possible choices of c. Thus, the asymptotics series of 
u>o is c-independent [0] and is given by 

5 Up to trivial rescalings of /q, ui and t as in (|6l|). 
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„ = ]T4 0) ^ (5 ^ 1)/4 (68) 
o 

where now there are even and odd contributions to the free energy as an 
expansion of integer powers of 1/N. It is also straighforward to rewrite the 
expression of the free energy 0, in the presence of a cubic interaction, 

4> = fo + u' (69) 

Let us finally remark that when c = 0, the mapping between uq and fo 
turns out to be an ordinary Miura transformation [15]. 
6. Concluding Remarks 

There is a good understanding of the features exhibited by the Hermitian 
ensemble in its several versions such as the one-matrix model, multimatrix 
extensions [|l6j |17| and the Hermitian matrix model in the presence of an 
external field (the Kontsevich model) |18" |, [19] , [p0[| . However, if we exchange 



the Hermitian ensemble with the real symmetric ensemble in the foremen- 
tioned models, our knowledge decreases drastically. Even in the simplest case 
of a real symmetric one-matrix model, our understanding is not that good. 
This makes it quite interesting to check the properties of the real symmet- 
ric matrices in comparison with the equivalent properties for the Hermitian 
matrices, as has been done in this letter for the case of pure gravity. 

Following this philosophy, it would be quite interesting to know whether 
the partition function of the real symmetric model is the r-function of a cer- 
tain hierarchy of integrable differential equations. The Miura transformation 
( |67D might be helpful for this purpose. Unfortunately, this is not the same 
Miura transformation relating the KdV hierarchy associated with (d 2 + /o) 
to its partner known as the mKdV hierarchy [[/J. 

There is a natural extension of the present work to check whether or not 
there exits a doubling of the ordinary differential equation flB"T|) when a more 



general non-even potential is considered ||21|| . The simplest case in which this 
feature should appear is with a potential like V(x) = g2X 2 + g^x 3 + (74X 4 , as 
it happens in the Hermitian model [TT|] fl2| |. 
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Appendix I 

Let us analyze the other possibility for the equation (|36|) . If a 2 k+i,2k 7^ 0, 
we end up with the following solution for the element a2fc+2,2fc+i, 

2 h 2k 2k + 1 

fl2fc+2,2fc+l — -"5 "5 {(V) 

pgz r k f)p2k+i 

But let us notice that this is precisely an element of the type a 2k>2 i~-i 
already determined in equation Q55p. Thus, comparing the equations ( |70|) 
and (|35|) we obtain now a recursion relation between r k and h 2k 

, 93 <f3 2 93 r k+l 2k + 1 

h 2k = — r- [-=-T + ] r k (71) 

z z n 2k+ i p 2k +i 

However, a 2k+ i t2k have not been determined yet. To do this, we have to use 

the remainig equation for j = 0. In this way, the roles of equations (Pi]) 

j = 1 and j = have been exchanged. 

It is convenient again to split indices in ( |i8|) j = in evens and odds. If 

i = 2k, as < R 2k -i, R2k >r= and < R 2k+ i, R 2k >= —r k , we have 

3 2k + 1 

h 2k = - [g3a2k+2,2k+l - -7, ](~ r k) (72) 

Z PP2k+l 

and this determines the last coefficient that we were left with being 

2 h 2k 2k + 1 

0-2fc+2,2fc+l — 3 1 a y' 6 

93P r k 93Pp2k+i 
Now, comparing the equations ([70|) an (f73|) we draw the conclusion that 
the possibility we are dealing with is only consistent iff the coupling constant 
takes the value 
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g 3 = -i (74) 

In addition we have still to impose the relation coming from equation (|18|) 



j = when the indices are odd. Unlike the case of the simple solution (|37j), 
now this relation takes an extremely long and cumbersome form due to the 
fact that most of the a-coefficients do not vanish. Nevertheless, as the cou- 
pling constant is absolutely fixed by (|T4]), the possibility (|70|), if consistent, 
would not be relevant as far as the contunuum limit is concerned. 
Appendix II 

To obtain the recursion relation (PO) solving the model, we start with 



equation fll8|) for j = 0. After some algebra involving the use of the obser- 
vations described in section 3, we arrive at 

hi = — 7^fi,i+2 < Ci + 2, Ri >R 

~ 2" u i,i+i < C i+ i,Ri >r ——Vij-i < Cj_i,i2j >r (75) 
Using ([16]) to invert the relation between Ri and Cj, we find 

hi = ~[v i,i+2di+2,i+l ~ Vi,i+l] < Ri+l, Ri >R 

— ^{^i,i+2[ — ( a i+2,i+l a i+l,i — a i+2,i)0"i,i-l + (0'i+2,i+l a i+l,i-l ~ Oj+2,i-l)] 

+ Vi,i+l( a i+l,t a i,i-l — aj + i ji _i) — fj^-i} < Ri-i, Ri >r (76) 

It is useful again to split the indices into even and odd. If % = 2k, as 
< R 2 k-i,R2k >r= and < R 2 k+i,R2k >r= it follows 

P r k+l /--v 

a2fc+2,2fc+l — 77^3 T {'<) 

* ft 2fc+l 

Replacing the value of 02^+2,2^+1 from (|35|) we obtain the desired result (fll|) 
for the recursion relation, 

ta* = h-j-^sT + 75- W (78) 

4 ^2fc+l ^P2fc+1 
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If i = 2k + 1, the analysis in principle is more cumbersome. Now we have 
< R 2 k+2,R2k+i >r= and < R 2 k,R2k+i >r= r k . Fortunately, using the 
equations (f28"|), fl3"2|) and (^) it easy to see that all the constants a's entering 
the equation (fFfl) are now zero, except for a i+ 2i-i = cJ2fc+3 2fc = — f^ir^ 1 - 
Then, after using (|19D, the equation takes the form 

/? 2 2 r fc+ i 2A; + 1 

Dividing ([FjJ) by p2fc+i = we achieve precisely the equation ( [78]) again. 
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